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O 

To observe CP-violating asymmetries through the interference of a weaker 
Q ■ amplitude with a stronger one in B° — > D^ix and B — > D^it decays, one 

must collect enough events that the intensity associated with the weaker 
amplitude would be statistically significant. We show that provided the 
weaker amplitude is measured separately in B^ — > D^^ir decays, the 
j> ■ time-integrated approach requires around 2.5 • 10 8 BB pairs for measure- 

0\ . ments of the weak phase sin(2/5 + 7) with an uncertainty of 0.05 or better. 

We also determine the optimal conditions for precise 2/3 + 7 measurements 
and discuss the possibilities for resolving a discrete ambiguity. 

PACS Categories: 13.25.Hw, 14.40.Nd, 14.65.Fy, 12.15.Hh 

I Introduction 

The phases of elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix describ- 
ing the weak charge- changing interactions of quarks are of fundamental importance. 
Together with magnitudes of the matrix elements and masses of the six quarks (u, c, t) 
and (d, s, b), these phases must be explained by any theory which extends our knowl- 
edge beyond the Standard Model (SM) of electroweak and strong interactions. 

Indirect information on CKM phases [TI 0] is now being supplemented by mea- 
surements of CP-violating asymmetries in B decays [§, f|, H which provide direct 
phase measurements. The weak phase (3 = Arg(— V* b V c d/V t * b V tc i) is determined by 
measurements of the rate asymmetry in decays such as B° — > J/ipKg, while a = 
Arg(— V t * b V t d/V* b V u d) will be determined by measurements in decays such as B — > nn 
and B — ► pn. Information on all charge modes will be needed to separate contributing 
amplitudes from one another 0. 

Information on 7 = Arg(— V* b V u d/V* b V c d) is more difficult to obtain. The decays 
B ± -> D ^, B ± -> U K ± , and B ± -> Dcp^, where D CP is a CP eigenstate, 
permit one to perform a triangle construction to extract the weak phase 7 0. The 
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interference of the Cabibbo-favored decay D° — > K~tt + and the doubly-Cabibbo- 
suppressed decay D° — > K + n~ introduces an important subtlety in this method [§]. 
Numerous determinations of 7 using nonstrange and strange B decays to tttt and Ktt 
are subject to questions associated with SU(3) flavor violation, electroweak penguin 
contributions, and rescattering 0. 

The Cabibbo-favored decays B° — > D^~ir + and B° — > D^ + ir~ and the corre- 
sponding doubly-Cabibbo-suppressed modes B — > D^~n + and B° — > D^ + K%~ 
can provide information on the weak phase 2(3 + 7 B, ffDl O, O. (One can substitute 
p ± or af for the charged pion.) These methods typically require measuring either a 
very small rate asymmetry (for the Cabibbo-favored modes) or a very small rate (for 
the Cabibbo-suppressed modes). It was therefore suggested recently JTB) that one 
instead measure 2(3 + 7 via the interference of a small amplitude with a larger one in 
decays of the form B — > V1V2, where, for example, Vi = -D* and V2 = p. The interfer- 
ence is to be detected through characteristic angular distributions in decay products 
of the vector mesons, and through time-dependent measurements. Refs. [[12] and [JTJ] 
contain some useful results regarding these distributions. 

In the present paper we analyze the possibilities of precise measurements of 2/3 + 7 
for the simplest case of B — > D^n decays. We find the optimal conditions for measur- 
ing 2(3 + 7. We also estimate the number of BB pairs needed for such measurements 
that will reduce the allowed range of 2(3 + 7 values to the currently achieved indirect 
bounds coming from measurements of other CKM parameters. 

A general feature of CP- violating asymmetries detected through the interference of 
a weaker amplitude with a stronger one is that one must be able to detect processes at 
the level of the absolute square of the weaker amplitude [15]. We find that this situation 
holds for B — > D^tc decays. One still has to be able to collect enough events such 
that the absolute square of the Cabibbo-suppressed amplitude would be detectable 
with good statistical significance. This translates to the need for several times 10 8 
produced BB pairs. (Ref. |1| cites a figure of 10 s pairs for a useful measurement 
of sin {2(3 + 7) using B — > ViV 2 decays.) In fact, our best determination makes use 
of a direct measurement of the weaker amplitude through a factorization relation 
between B° — > D^ + ix~ and B + — > D^ + n° |H|| . For both pseudoscalar and vector 
D mesons in the final state, we employ different models to anticipate the size of the 
weaker amplitude. However, direct measurements of the rates for B + — ► D + tt° and 
B + — > D* + 7T° will eventually give us these amplitudes directly. 

In Section II we introduce our notation and predictions for decay rates of neutral 
B mesons in the framework of the time-integrated approach. We shall quote results 
for B — > D*tt decays because of advantages in D* detection, recognizing that many 
are also valid for B — > Dir. Decay rates as functions of a minimum vertex separation 
(expressed in terms of proper time) are of particular interest to us in Section III as 
we try to find the optimal conditions for measuring the weak phase 2/3 + 7 with high 
precision. In Section IV we circumvent the problem of measuring the small weaker- 
to-stronger amplitude ratio R by making a foray into charged B meson decays, using 
the process B+ -> ZJM+tt . Estimates of the minimum number of BB pairs required 
for precise measurements of 2(3 + 7 are obtained in Section V. These are convoluted 
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with a finite time-resolution function and realistic mistagging probabilities in Section 
VI. In Section VII we discuss a possibility of partial resolution of an 8-fold discrete 
ambiguity by separating 2/3 + 7 an d the strong phase S between Cabibbo-allowed and 
Cabibbo-suppressed modes. We summarize our results in Section VIII. 



II Notation and predictions 

The "right-sign" decays B° — > D*~n + and 5° — > D* + tt~ are governed by the Cabibbo- 
favored combination of CKM matrix elements V* b V u d or charge-conjugate, while the 
"wrong-sign" decays B° — > D*~n + and B° — > D* + n~ are governed by the doubly- 
Cabibbo-suppressed combination V* d V u b or charge-conjugate. We denote / = D*~n + 
and / = D* + 7i~ . Then from 

(f\B°) = A 1 e i + 1 e iSl , (f\B°) = A 2 e l<t,2 e l&2 (1) 

it follows that 

(f\B°) = Axe'^e 151 , (f\B°) = A 2 e-^ 2 e i52 , (2) 

where the weak phases (pi change sign under CP conjugation, while the strong phases 
Si do not. The amplitudes are in the ratio 



„ A 2 
Ax 



v; d v ub 



v; b v ud 



- A 2 (p -irj)\r ~ 0.02r , (3) 



where A ~ 0.22, p, and rj are parameters |L6[ which describe CKM matrix elements 



and r = 0(1) describes a ratio of decay constants and form factors. The weak phase 
difference is 

0i-0 2 = Argf^) =7T + 7 . (4) 

\ V cd V ubJ 

We write the time-dependent decay amplitudes in terms of the functions 

U(t) = e- imt e- rt/2 cos(Amt/2) , /_(*) = e~ imt e - rt/2 t sin(Amt/2) , (5) 

where m — (rriL + mu) /2 is the average of the two mass eigenvalues, Am = m# — iul 
is their difference, T = (Ti + T 2 )/2 is the average decay rate of the eigenstates, and 
we neglect Ar = T H — T L . Then 

(f\B°(t)) = f + (t)(f\B°) + p.(t)(f\B°) 

= e- mt e- rt / 2 (A 1 e i ^e^cos— + tlA 2 e^ 2 e l52 sm—) , (6) 

V 2 p 2 J 

(f\B°(t)) = / + (*)</|B°> + |/_(*)</|B°> 

= e- im V r */ 2 fA ie -^e^cos^ + z^ 2 e-^e i52 sin^) .(7) 

V 2 q 2 J 



If B°-B mixing is described primarily by standard model loop contributions domi- 
nated by intermediate tt pairs, we have q/p = e~ 2 ' il3 , and 



(f\B°(t)) 



A 2 
2 



1 + R 2 + (1 - R 2 ) cos Amt 



- 2R sin (2(3 + 7-5) sin Amt] 



(fmt)) 



A 2 

2k p -rt 
2 



1 + R 2 + (1 - R 2 ) cos Amt 
+ 2R sin (2(3 + 7 + 5) sin Amt] , 



(9) 



where 5 = 5 2 — 5i. 

Retracing the above steps for the "wrong-sign" decays B° — > D* + n~ and B 
D*~n + , we find 



(f\B°(t)) 



(fmt)) 



-Yt 



1 + R 2 - (1 - R 2 ) cos Amt 



- 2R sin(2/3 + 7 + 5) sin Amt] 



2 



1 + i? 2 - (1 - R 2 ) cos Amt 
+ 2i? sin(2/3 + 7 - 5) sin Amt] . 



(10) 



Let us now consider the production of a B°B° pair in e + e — > T(4S) — > B°B^, so 
that the pair is in a state of negative charge-conjugation eigenvalue. Assume that 
we "tag" the initial production of a B°(p) with a B°(—p), and the initial production 
of a B°(p) with a 73 (— p). Then if we define the proper decay time of the state / 
with center-of-mass direction p as t/, that of the tagging state with direction — p as 
t t , and t' = t f - t t , T = t f + t t , we find @, p], 



Ki^-^zr^)!^)! 2 
|(^(-p),d*Vf(?0|¥->| 2 



-IT 



1 + R 2 ± (1 - R 2 ) cos Amt' 



2R sin(2/3 + 7 + 5) sin Amt'] 



-FT i 



l + i? 2 ±(l -R 2 ) cos Amt' 



(12) 



+ 2R sin(2/5 + 7 + 5) sin Amt'] . 
One can express the time-integrated decay rates as 

/ dt f dt t l(B\-p),D^(p)\m_)\ 2 oc / dt'e^A^f) + B^(t')} , 
f 00 dt f l°° dt t \(B°(-p),D* ± ^(p)\^^\ 2 oc [°° dt'e- r ^[A±(t')-B ± (t')} (13) 

JO JO J-oo 



where 

A±(f) = (1 + R 2 ) ± (1 - R 2 ) cos Amt' , 
B± (f) = -2R sin(2/3 + 7 + 5) sin Amt' 
are even and odd functions of t', respectively. 



(14) 
(15) 
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Now we introduce notation for measurable decay numbers. The number of B° 
D*~n + decays with vertex separation t' > is 



Nl<x dt'e~ m [A + {t') + B(t% (16) 
Jo 



while those with t' < is 



/0 p oo 

dt'e- m [A + {f) + B_(t')] = / dt'e- r W[A + (t') - B_(t% (17) 
-oo JO 

Here the superscript "r" denotes right-sign decays. The corresponding expressions 
N™ for the wrong-sign (superscript "w") decays B° — > D* + tt~ with £' > and iV^ 
for 5° -> D*+7T- with f < are 

AToc/ rft'e- r| * l [A_(t')±5+(t')]. (18) 

Similar expressions for 5° decays are 

iVloc / ^^(OTfl+^l. (19) 
Jo 

A^oc/ dt'e- m [A(t')^B(t')}. (20) 
Jo 

Note that the following 4 linear relations among the 8 decay numbers 



N r, W + ^ = ^ + N r,w^ (21) 
^ _ N r, W = J^r _ J^r ^ ^ 

limit the number of independent quantities to 4. In principle, that allows one to 
forgo measurements of B decay numbers. However, that method would lead to 
larger uncertainties in determination of 2f3 + 7 and we shall not use it. 

We shall investigate the dependence of the time- integrated rates on a minimum 
vertex separation to- The aim of the calculation is to find the optimal conditions 
for measuring sin(2/5 + 7). Fig. 1 shows that indirect bounds on that weak phase 
coming from measurements of other CKM parameters p, [Tj| limit the expected value 
of sin(2/5 + 7) to the region between 0.89 and 1. To get in the same ballpark with 
the indirect bounds we will calculate the number of BB pairs required to determine 
sin(2/5 + 7) with an uncertainty of 0.05. This is the main goal of the paper. 

Ill Decays with vertex separation greater than to 



If one only takes into account decays with vertex separation greater than to, Eqs. ([16- 
0[) become 



POO 

N r ± {t )oc dt' e- m [A + {t')±B_(t% (23) 
J t 

AT- (t ) oc / dt'e-^ l[4_(0 ± B + (t% (24) 
Jtn 
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Figure 1: Contours of sin(2/3 + 7) (thin curves with values to the right) in (p,rj) plane. 
Thick lines denote current limits on CKM matrix parameters 0, |18| . Solid circles 
denote limits on IV^/V^I from charmless b decays, dashed circles denote limits on 
Vtd from B° — B mixing, and the dotted circle denotes the lower limit on iV^/Vy 
from the lower limit on B s — B s mixing. Dot-dashed hyperbolae come from limits on 
CP- violating K° — K mixing (the parameter e). Two solid rays denote the recent 
world average ±lcr limits sin(2/3) = 0.79 ± 0.10 from neutral B meson decays. The 
allowed range is shaded gray. 
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roo 

NUto)<x / dt'e-r^iA+WTB+it')}, 

Jt 

A^(t ) oc / dt'e' m [A4t') =F £_(*')]■ 



(25) 
(26) 



There are several ways to combine these decay numbers together into algebraic sums. 
Some of the resulting combinations may include one of the following expressions: 
A+(t') + A_(t'), A+(t') - B+(t') + B_(t'), or B + {t') - £_(*')• Composing 

ratios of these algebraic sums (see fi(t ), /^(to) and fzito) below), we can extract the 
parameters R, sin(2/3 + 7) cos 5 and cos(2/5 + 7) sin 5: 



R 



SC = sin(2/5 + 7) cos 5 
CS = cos(2/? + 7 ) sin 5 



a(t ) - A (t ) 
\a{t Q ) + f l (t y 

1 + R 2 



2b(t ) R 
1 + R 2 



h(t ) 



where 



26(* ) R 

[•CO / I 

a(t ) = Te rto / dt'e~ rt cos(Amt') = cos^rto + A), 

/•oo I 

6(*o) = re r *° / dt'e~ rt sin(Amt') = —j= sin(x d n + A), 
J* Jl + x 2 



and 



fi(to) 
/ 2 (to) 
/b(«o) 



A = arctanxd, 


x d = Am/r, 




(iv; + jvi + ni + nz) 


- (JV™ + J\T + A7f 


+ N?) 










- (Nl + N™ + W 


+ NlO) 




N 




(iV; + iV™ + Nl + ~N») 


- (NL + N™ + NL 




N 



with 



(27) 

(28) 
(29) 

(30) 

(31) 
(32) 

(33) 
(34) 
(35) 

(36) 



iV = Nl + NL + Nl + NL + N™ + N™ + N% + N™. 

We have suppressed (to) after the decay numbers in the last four formulae. 

It has been noted in [13| that R is too small to be determined by this method. 



Indeed, calculations show that the smallest uncertainty in R is achieved at to — and 
is equal to 

1 4 (2 + x* d ) 1 _ no 



a(R) 



■ 0.03, 



\| 16R 2 e(B r + B W )N B 

with e being the tagging efficiency We take e to be 0.684±0.007 M. B r , the branching 
ratio of the "right-sign" decays B° -> D*-<ir + , equals (2.76 ± 0.21) x 10~ 3 |9|. One 
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r (D*ir) 


r (£>7r) 


Light Front Model 


[22j 0.81 


0.72 


BSW II Model 


[ 


?! 


| 1.33 


1.11 


NS Model 


24 




0.88 


0.72 


LCSR Model 


25 


1.01 


0.87 


MS Model 


1 




1 


0.92 


0.82 



Table I: The ratio r evaluated in various models. 

can show that for Xd = 0.756 ± 0.012 |2C| the branching ratio of "wrong-sign" decays 
is B w = kB r » 0.61 ■ lO" 3 , with k « x 2 J{2 + xg) « 0.22. 

The error a(R) = 0.03 is bigger than the approximate R value itself [Eq. (|3|)]. 
Thus, one has to search for another method of measuring R. 



IV Ratio of amplitudes 



The main reason one cannot get R directly from the ratio of B° — > D^ + ix and 
B° — > /} ( *)~7r + decay rates is that the large 5° — > B° mixing amplitude in the former 
overwhelms the smaller direct tree contribution. One can circumvent this obstacle 
by considering decays of charged B mesons, e.g. B ± — > D^^tt , as suggested in [|TT|]. 
The tree amplitude is dominant in these decays and is proportional to A\/2. Thus, 
the ratio of B ± — > Z)'*^/ and B° — ► D^*^n + decay rates can be used to provide a 
simple way to estimate R. 

The B + -> D ( * )+ 7T° decay rate can be estimated by assuming factorization: 



M = ^fv: b V cd (n(p 



-q)\h,u\B(p))(D^\V,\0) 



(3? 



Using the standard parameterization | 21| , one obtains the ratio r defined in Eq. 



r{Dir) 



Id* F{ 



Btt, 



m 



D* 



f D (m 2 B - ml) F£ 



Btt, 



m 



D) 



fAm 2 B-ml)F BD (ml] 



(39) 



In Table I, we give the values of r for B° — > D < ^*' ) 7r decays in several models. In all 
cases, the models predict that r is close to unity, i.e. R ~ 0.02. 

The error on R can be estimated using the method described in the beginning of 
this Section. Suppose that the number of detected B + — ► D^+ir decays is N 2 out 
of iV tagged I? +, s, while the number of detected B° — > Z}(*)~7r + decays is N\ out of 
the same number N of tagged B 0, s. Then, assuming equal charged and neutral B 
production, the value 



R 



!— 



iVi 



(40) 



S 



has an uncertainty 



a(R) 



1 B 



(41) 



Taking into account B — > tt° and B° — > D^*^ + 7r decays increases statistics 
by a factor of 2, leading to cr(-R) = l/^/2N r , where N r = 2N\ is the number of 
B° — > D(*^7r + decays plus the number B° — > /}W + 7r~ decays. To make connection 
with the total number Nb of produced BB pairs, note that the number of tagged 
events is N r = eB r Ns- Thus, 



a(R) 



V2eB r N B 



(42) 



For 10 8 produced BB pairs (j(R) = 0.17 ■ 10~ 2 , i.e. less than 10% of its value. Thus, 
measurements of B + — > D^^n decay rates provide the ratio of amplitudes with 
a high precision. This information may be used in the time-integrated approach 
discussed in the previous Section. Now we can go a step further and estimate the 
uncertainty in determination of sin(2/? + 7) cos 5 and cos(2/3 + 7) sin S. 

In the following analysis, we will take r = 1 (corresponding to R ~ 0.02) and use 
Eq. ([4"2"D to estimate the error on the ratio R. 



V Uncertainties in sin(2/3+7) cos 5 and cos(2/3+7) sin 5 
with perfect time resolution and no mistagging 

The uncertainties in the ratios f 2 and f 3 [see Eqs. ( p4|) and fl35|)1 are 

1 



iV (to) ^ + B-)iV fl \ 6(1 + fc) BtJVb 



(43) 



Eqs. (|28|) and fl2"9|) allow an estimate of the values of /2 and fa: f 2 ~ 2Rb(t ) SC, 
« 2Rb(to) CS. Now that Eq. ([|2|) provides the error in i?, we can calculate the 
uncertainties in S'C and CS 1 : 



<r(SC) 



1 



R 4 



-a 



] (R) + 



< 



1 



2(1 + k) 6 2 (t ) + e r *° 



R 2 -2b{t )R\ e(l + k)B r N B 



(44) 



2b(t 



i? 4 



i? 2 26(t ) i? 26(t ) \ e(l + k) B r N B 



(45) 

Finally, one can calculate the number of BB pairs needed to get any particular 
precision a (SC): 

2(1 + k)b 2 (t ) + e Tt ° 



4e(l + k) b 2 (t )R 2 B r a 2 (SC)' 



(46) 
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Figure 2: Number of produced BB pairs needed to achieve an uncertainty of 0.1 in 
measurements of sin(2/3 + 7) cos 5 (solid line) and cos(2/5 + 7) sin<5 (dashed line) vs. 
minimum vertex separation t . Perfect time resolution is assumed. 



or a {CS): 



D Ft 



Nb ~ 4e(l + fc) b 2 (t )R 2 B r a 2 (CS)' ^' 

As seen from the figure, these two quantities have the same minimum location because 
they only differ by a constant independent of to- Here we have assumed that fz{to) 
is proportional to cos(2/9 + 7) sin 5, which is expected to be small, and that SC is 
close to I. However, the neglected SC and CS dependence can be readily put back 
if necessary, and one finds that the position of the minima would remain the same, 
independent of the values of SC and CS, for both curves. 

Fig. 2 shows the Nb dependence on to according to the above two equations. The 
curves were calculated under the assumption that one needs to get <xo = 0.1. We found 
out that this precision level is sufficient to determine sin (2/5 + 7) with an uncertainty 
of 0.05 (Section VII). The optimal conditions for both measurements are achieved 
if one only takes into account decays with vertex separation greater than ~ 0.45/r. 
That one needs fewer BB pairs to reach the same precision for cos(2/3 + 7) sin 5 as 
indicated in Fig. 2 reflects our previous assumption of small cos(2/3 + 7) sin 5. Thus, 
the minimum uncertainties one can obtain if N B BB pairs are available are 



a mm {SC) ~ O.I1 



'1.62 • 10 8 



(48) 
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/O 96 • 10 8 

a min (CS)^ O.iy ' — . (49) 

Now we shall check how these formulae change if we take into account finite time 
resolution and realistic mistagging probabilities. 



VI Finite time resolution; mistagging 

Measurements of the decay numbers are smeared by finite resolution of vertex sep- 
aration. For simplicity we shall assume a single Gaussian resolution function. The 
observed decay numbers are given by Eqs. (|23) — (p6|) convoluted with the resolution 
function 

/+°° I / , , 

d^^^e-^ 2 ' 2 " 2 . (50) 
-oo v Zna 



For example, 

N±fo)=N R(t )®N5 : (ji), (51) 

and similar convoluted relations for N±(t ), N±(t ), and N±(t ). Here N is a nor- 
malization factor and a is the resolution of time separation between vertices. For the 
BaBar detector the average resolution of space separation between vertices is 180 /im 
[fH while the average separation is (3^cr B o = 260 /xm, implying aT = 180/260 = 0.69. 

The algebraic sums of decay numbers that enter Eqs. (Bl) and (RSI) have to be 
modified correspondingly. For example, Eq. ( |36|) becomes 

n M = 2 ^£y^(^jp dt > e -m [A+(tl) + A _ V)] 

Nn f+°° ~ ,r, r 1 

= 2-M dt e" 1 * 1 \A + (t/T) + A_(t/r) 

1 J —oo 



= 4^ (l + |)(J 1 + J 2 ), (52) 

Similarly, one obtains 

(NL+N™ + N$+N%) - {N^ + N^+W + T^) = 8^ J RJ 3 5C, (53) 

(Nl + + NI + W!) - (N r _ + +W + Nf) = 8^ RJ 3 CS. (54) 

In the above equations, $(x) = (2/-y/7r) J X e~ z2 <i2 is the error function and 

/■+oo 

Ji = / dt = 2, (55) 



J 2 = / + °°^ *(^ £ )e HAl , (56) 



J 3 = jTj'd/i ^(t^-] e -^ s inx d jl, (57) 
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The last two integrals have been numerically evaluated for different values of to i n 
the range from to 1.5/T. Now SC and CS can be rewritten in terms of the ratios 
f 2 and f 3 as 

SC,CS=^±^^^f 2 M. (58) 

Next, we will take into account the mistagging factor. Mistagging refers to the 
cases where a decay (B° — > tag, B° — > D*~n + ) was incorrectly identified as (B° — > tag, 
B° — > D*~7i + ), and vice versa. Thus, one sees that decays labelled as B° — > D*~tt + 
("right-sign" decays) actually contain some 5° — > D*~n + ("wrong-sign" decays) 
events. As a result, experimental measurements only provide decay numbers smeared 
by the mistagging effect. For instance, the numbers of apparent right-sign events are 



w)N r ± (t )+wN%(t ) 



(59) 



where w is the mistagging probability. For the BaBar detector the tagging efficiency is 
e = e « — 0.684 ± 0.007 while the effective tagging efficiency is Q = e i(l ~ 2wi) 2 = 
0.261 ± 0.012 Q. For our purposes we will simplify calculations by assuming the 
single tagging option with e = 0.684 and Q = e(l — 2w) 2 = 0.261. Thus, the effective 
mistagging probability is w = 0.191. 

Note that the sum of all smeared decay numbers is still equal to N, the sum of all 
physical decay numbers. One can show that the ratios f' 2 and f' z composed of smeared 
decay numbers are related to f 2 and f% by f 23 = (1 — 2w)f 2;3 . Thus, experimental 
measurements of smeared decay numbers allow the direct calculations of SC and CS: 



sees 



J 1 + J 2 1 + B 2 



l-2w 2 J, 



R 



i 2,3(^0 



(60) 



Assuming that experimental uncertainties are <?[N±(to)] = JN±(to), ff[N±(to)] = 
N±'(t ), etc., we can estimate the uncertainties in f 2 and f 3 measurements to be 



1 



fW(^ \e{l+k)BrN B - 
The uncertainties in SC and CS measurements are 
1 Ji + J 2 1 



(61) 



cr(SC) < 



l-2w 2 J 3 R\ 



2(1 + k) [J 3 /( Ji + J 2 )] 2 (1 - 2wf + e r 'o 



e(l + k) B r N B 



a{CS) 



1 J1 + J2 1 



(62) 



(63) 



l-2w 2J 3 R\\e(l+k)B r N B ' 

We assumed a small CS in deriving the second equation. Finally, one can calculate 
the number of BB pairs needed to get any particular precision <jq(SC): 



N. 



B 



2(1 - 2w) 2 (l + fc) Jj + e lt0 (Ji + J 2 ) 
4e(l - 2w) 2 (l + k) Jf R 2 B r a 2 (SC) 



(64) 
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Figure 3: Number of produced BB pairs needed to achieve the uncertainty of 0.1 in 
measurements of sin(2/3 + 7) cos 5 (solid line) and cos(2/? + 7) sin 5 (dashed line) vs. 
minimum vertex separation t . 



or a (CS): 

N e rt »(J 1 + J 2 ) 2 

B ~ 4e(l - 2w) 2 (l + jfe) J| IP&oKCS) ' 1 J 

As in the previous Section, the position of the minima is the same for both curves 
and is independent of the values of SC and CS. 

Fig. 3 shows the N B dependence on t . The curves were calculated under the 
assumption that one needs to get cq = 0.1. The optimal conditions for measurements 
are achieved if one only takes into account decays with vertex separation greater than 
0.44/r. Then 

, QrA ni / 5.06 ■ 10 8 
(j m in{SC) ~ 0.1W — — — , (66) 



/4 40 • 10 s 

a min {CS) ~ O.lW Nb ■ (67) 



If BaBar is able to improve its performance to the level quoted in [27], i.e. a(Az) = 
110 /im, e = 0.767 and Q = 0.279, then the required minimum number of BB pairs 
reduces by a factor of ~ 1.4 for both SC and CS measurements. Besides, the position 
of the minima is shifted to a slightly larger value (t Q ~ 0.53/T) of vertex separation. 
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Figure 4: Contours of s_, s + , and their uncertainties cx(s_) and cx(s + ) in the (SC,CS) 
plane. Only the first quadrant of the plane is shown. The plots in other quadrants 
are symmetric to those in the first one since the plotted quantities only depend on 
the absolute values of SC and CS. The blank triangle above the line SC + CS — 1 
denotes the forbidden region on the plane: SC + CS = sin (2(3 + 7 + 5) should always 
be smaller or equal to 1. When CS = cr(s-) achieves its smallest values: 0.1 and 
0.07 for plots (c) and (e) respectively. 
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VII Extraction of sin(2/3 + 7) and cos 5 



If one measures sin(2/3 + 7) cos 5 and cos(2/5 + 7) sin 5 values to be SC and CS, then 
trigonometry dictates the following values for sin 2 (2(3 + 7) and cos 2 5: 

sin 2 (2(3 + 7), cos 2 6 = sr% = ^(l + SC 2 - CS 2 ± y/\(l,SC 2 ,CS 2 )^) , (68) 

where X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2yz — 2zx. When one root corresponds to 
sin 2 {2(3 + 7), the other corresponds to cos 2 5. There is an ambiguity: Which is which? 
One cannot resolve it without making additional assumptions. 

If the value of sin (2(3 + 7) is assumed to be in agreement with the indirect bounds 
(Fig. 1) then it should be larger then 0.89. However, indications that cos 5 is large, 
too [E3], do not allow an easy distinction between the two quantities. Figs. 4(a,b) 
show the contours of s_ and s + values in the (SC,CS) plane. There is a big region 
on the plane where s_ < 0.89 while s + > 0.89. If the measured values of SC and 
CS fall inside this region then sin (2(3 + 7) could only be associated with s + and 
the ambiguity might be resolved. The possibility of resolution also depends on the 
uncertainties cr(s_) and a(s + ). Those are calculated from Eq. ( |B8D with the help of 
Eqs. (|66D-(|67|). The contours of these uncertainties are shown in Figs. 4(c,d,e,f) for 
two different numbers of BB pairs. For example, if the number of produced BB pairs 
is 5.06 ■ 10 8 and (SC, CS) = (0.75, 0.15), then we can calculate s+ = 0.97 ± 0.04 and 
s_ = 0.77±0.12. In this case, s_ does not take values that are larger than 0.89 within 
the la level, and the solution favored for consistency with Fig. 1 is sin 2 (2(3 + 7) = s 2 + , 
cos 2 5 = s 2 . The 4-fold ambiguity in 2(3 + 7 remains but reduces to a 2- fold one 
when we take into account that only positive values of sin(2/? + 7) are consistent with 
indirect bounds. One can see from Fig. 1 that if 2(3 + 7 < tt/2 then sin(2/5 + 7) 
should be larger than 0.97. This fact might completely resolve the ambiguity in favor 
of it/2 < 2(3 + 7 < 7r if values larger than 0.97 are measured to be inconsistent with 
sin(2/3 + 7) within the la level. 

Of course, one cannot exclude the possibility that sin(2/5 + 7) is inconsistent with 
indirect bounds and is substantially smaller than 0.89 while | cos<5| is close to unity. In 
that case, one would make a wrong assignment of s + and s_ to sin(2/3 + 7) and cos S, 
respectively. Therefore, it is preferable to make other measurements of sin(2/5 + 7) 
in decays like B — > D^p or B — ► D^ai where strong phase might differ from 5 in 
B — > D^n decays. 

It is also worth noting that for the overwhelming part of the region where s+ > 
0.89, the uncertainty in sin(2/3 + 7) is at most 0.05 [cf. Figs. 4(b) and 4(d)]. Thus, 
for many values of SC and CS this method allows a very precise determination of 
sin(2/5 + 7) and a good measurement of the strong phase S. 

Besides, the method can be used to detect deviations from the Standard Model. 
If the measured values of sin(2/3 + 7) cos 5 and cos(2/9 + 7) sin 5 fall into the upper 
left corner of the (SC,CS) plane, then both s_ and s + would be inconsistent with 
the 0.89 — 1.0 range expected from the unitarity of the CKM matrix. 
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VIII Conclusions 



This paper has explored the optimal conditions for measurements of weak phase angle 
2(3 + 7 and strong phase 5 between Cabibbo-allowed and doubly- Cabibbo-suppressed 
amplitudes in B — > D^ti decays. We have found that in the time-integrated approach 
it is advantageous to only consider events with vertex separation greater than t 
which is equal to 0.44/r for the BaBar detection parameters. The loss in statistics is 
outweighed by an increase in the integrated asymmetry. 

Fig. 3 shows that production of approximately 5- 10 8 BB pairs is needed to reduce 
the uncertainty in determination of sin(2/3 + 7) cos 5 to 0.1 in B — > D*ir decays. A 
smaller error on cos(2/3 + 7) sin 5 will be achieved at the same time if its value is small. 
B — > Dtt decays have the advantage of a slightly higher branching ratio but a setback 
in D meson detection. The combination of both types of decays might reduce the 
number of needed BB pairs to 2.5 • 10 8 , an amount within the reach of both BaBar 
and BELLE in the next few years. A time-dependent analysis fl7| does not lead to 
any improvement with respect to this figure. 

If the strong phase 5 is not very close to or tt, the ambiguity between sin(2/3 + 7) 
and cos S can be resolved. This method allows sin(2/3 + 7) to be determined with a 
precision of 0.05 or better. 
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